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ABSTRACT 


FAMSPJUB  is  a  digital  computer  program,  uritten  in  single  pre¬ 
cision  FORTRAN  IV/  for  the  Philco-Ford/GE  265  time-sharing  system, 
to  determine  the  eigenf requencies  and  eigenmodes  of  vibration  in 
pure  bending  of  loaded,  piece-uise  continuous  single-span  beams 
supported  in  any  one  of  the  ten  possible  combinations  of  the  fol¬ 
lowing  end  conditions :  free,  guided,  hinged  and  clamped.  Continu¬ 
ously  variable  beam  properties  may  be  accommodated  by  appropriate 
subdivision  of  the  span.  The  distributed  mass  of  the  beam  itself 
is  taken  into  account  by  means  of  field  transfer  matrices ;  this  re¬ 
presents  an  advance  over  the  usual  lumped-parameter  approach.  To 
preserve  numerical  accuracy  at  the  higher  mode  values,  a  modified 
transfer  matrix  method  is  employed.  Uerification  solutions  have 
been  obtained  on  the  IBK  360-50/65  using  the  double ,  precision 
facility. 
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FOREWORD 


This  report  documents  results  obtained  from  project  mork 
on  task  title  "In-House  Initiated  Lab  R  &  D"  (DA  Proj.  IMo. 
1TD61101A91A/Task  Area  DQ/bJork  Unit  IMo.  2D5EH)  authorized  by 
Dr.  E.M.  Patrick,  USATACOM  Chief  Scientist/Technical  Director 
of  Laboratories.  The  project  is  intended  to  strengthen  in-house 
research  capabilities  and  provide  a  methodology  for  dealing 
uith  practical  problems.  The  author  extends  his  sincerest 
appreciation  to  Dr.  Patrick  for  the  opportunity  to  perform 
this  uork  and  to  his  senior  colleague,  Mr.  Douglas  J.  Hack- 
enbruch  for  the  encouragement  and  assistance  he  has  provided 
during  all  phases  of  this  project. 


INTRODUCTION 


This  work  uas  initiated  after  a  literature  search  revealed 
that  no  generally  applicable  computerized  solution  for  beam  vibra¬ 
tion  eigenvalues  had  yet  been  developed  uhich  took  account  of  the 
continuity  of  beam  properties.  The  analytical  development  con¬ 
tained  herein  has  been  applied  to  piece-mise  continuous  beams  tak¬ 
ing  account  of  bending  only.  Extensions  uhich  mould  include  con¬ 
tinuous  variation  of  properties,  the  effects  of  shear  deflection 
and  rotary  inertia,  intermediate  support  conditions  and  response 
to  farcing  functions  are  quite  feasible. 

Conventional  numerical  methods  of  determining  natural  fre¬ 
quencies  of  transverse  beam  vibration  generally  employ  variations 
of  lumped-parameter  techniques.  To  account  in  a  more  refined  may 
for  the  distributed  mass  of  the  elastic  system,  field  transfer 
matrices,  uhich  mill  be  developed  in  the  follouing  section,  have 
been  incorporated  into  FAMSWUB  and  combined  uith  point  transfer 
matrices  to  accommodate  discrete  loading.  In  addition  to  discrete 
loading,  step-uise  variations  can  be  accommodated  in:  uniformly- 
distributed  loading,  beam  material  density,  elastic  modulus  and 
section  modulus. 

The  source  language  employed  should  provide  greater  time¬ 
sharing  system  compatibility  insofar  as  no  use  is  made  of  sub¬ 
scripted  subscripts,-  mixed  modes,  variable  name  lengths  exceeding 
six  characters  nor  of  the  double  precision  facility.  To  accommo¬ 
date  users  .  of  time-sharing  systems,  program  length  limitations 
have  been  observed  by  subdividing  FAMSIMUB  into  5  saved  items  deal¬ 
ing  uith  the  conditioning  of  input  data,  eigenfrequency  calcula¬ 
tion  (3  parts)  and  the  computation  of  normal  modes.  This  has  been 
accomplished  through  the  employment  of  intermediate  permanent  file 
storage  and  the  $USE  monitor  line. 

The  considerable  numerical  magnitude  difficulties  encountered 
in  evaluating  2x2  frequency  determinants,  uhich  involve  taking 
the  difference’betueen  quantities  uhich  are  several  orders  of  mag¬ 
nitude  larger  than  the  difference  itself,  are  sidestepped  by  an 
adaptation  of  a  method  developed  by  Pestel  &  Leckie  (1),  pp.  2B4- 
213. 

FAMSIMUB  commences  uith  the  first  non-zero  frequency  and  pro¬ 
ceeds  continuously  to  the  modeCs)  of  interest..  Houever,  a  start 
may  occur  at  any  frequency  to  investigate  a  resonance  condition. 
If  the  eigenmodes  are  also  desired,  a  printout  can  be  obtained 
in  four-column  format  of  deflection  amplitude,  maximum  bending 
stress  and  average  transverse  shear  stress,  together  uith  the 
corresponding  distance  along  the  beam  (Figure  27). 
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OBJECT 


Develop  a  uorking  method  for  computing  the  eigenvalues 
and  eigenmodes  of  transverse  beam  vibration,  taking  account  of 
the  distributed  beam  mass  in  preference  to  employing  a  lumped- 
parameter  approach. 
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Forces  and  Deflections  for  Beam  Element 
Figure  1 


Field  and  Point  Transfer  Matrices 

Consider  an  initially  straight,  piece-uJise  continuous  beam 
performing  free  transverse  vibrations.  Figure  1  shows  an  isolated 
beam'  element  of  length^  between  points  (i-l)  and  (i)  ,  having 
continuous  elastic  properties  and  a  constant  rate  of  inertia 
loading — ^  is  the  non-reactive  loading  per  unit  length  due 
to  beam  .self-weight  and  continuous  loading.  The  origin  of  the 
coordinate  plane,  in  which  all  vibration  is  assumed  to  occur, 


is  to  be  considered  as  coincident  with  the  left-^hand  support. 
The  sign  convention  is  also  revealed  by  Figure  1,  which  shows 
positive,  instantaneous  values  of  the  state  vector  components; 
deflection  w  ,  slope (in  radians),  bending  moment  M  and  shear 
force  V  .  The  uniform  bending  stiffness  for  this  length  is  Ejlj  , 
the  product  of  the  elastic  modulus  and  the  second  moment  of  area 
of  the  cross-section. 

From  simple  bending  theory,  the  fourth  partial  derivative 
of  w(x,t)  with  respect  to  x  is  proportional  to  the  rate  of  loading, 

thus  /  » 

El  a.  wM  =-a.w  Cl) 

Qx4  9 


For  a  solution  which  permits  eigenf requency  determination  we  write 
w  as  the  product  of  two  single-variable  functions 


w(x,t)»  W(x)sinwt  (2) 

whereW(x)  is  the  dynamic  amplitude  and  u  is  the  simple  harmonic 
circular  frequency;  substituting  into  (1) 


d^W 
dx 


which  becomes 


(3) 


d^W 

dx^ 


-m^W*5  0 


(4) 


(5) 
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Reuriting  (4)  in  LaPlace  transform  notation 


uhere 

solving 

ujhere 

Making 
f  ornis : 

W(x)  = 


W'(x)  = 


W"(x)  = 


W  ( p )  -  p*  Wq-  p’ wi-  pW  j;  -  W  o  -  W(  p )  =  0  ( 6 ) 

'W(p)=j^”-'“f(x)dx 

w(p)=  ^^[p’Wo+p’Wi+pW2+W'"j  (7) 

Wo=wf.,;Wj='J't,;  EiliWS^Mf^;  E;I^W^  =  Vil, 

the  appropriate  substitutions  and  applying  inverse  trans- 


W 


L_  cosmx +  coshtT>x sinmx  +  sinhmx 


2m 

j.  coshmx-cosmx  r^r  j.  sinhmx-sinmx ,,r 

2Ejljm2  2Eilim3 


(8) 


.t,l  _  m(sinhmx-sinmx)4A/R  4.  cos mx+ cosh mx  ^t/R 

„  ''•-I 


i.  sinmx  +  sinhmx  rrR  x  coshmx-cos mx  ^r 

x,-  .  ^i-i  X.-  .  9  ^  ^ 


2Eiljm 


2Eiljm2 


L  _  Ejljm^fcosh mx-cosmx) ^r  +  Ejljm(sinhmx-sinmx)  m^.R 


Mr  = 


i-i 


cosmx+coshmx  ..R  +  sinmx  +  sinhmx  ,,r 

- M:  ,  ■  Vj 


^I-I 

(10) 


2m 
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Transfer  Matrix  for  an  Elastic  Field 


yL  _  Eilim^(sin  mx  +  sinhmx)  yyR  +  Eilim^(coshmx-cosmx) \|/’R 


W"'(x) 


a.  m(sinhmx-sinmx)  4.  cos  mx  +  coshmx  .,r 

Z  '''‘i-i  o  ''i-i  ^ 


Hence  the  relation  of  the  state  vectors  Zj_^  and 


on  either  side  of  an  elastic  field  of  length  x  can  be  stated 
as 


Z- =  F  -Z- 

I  I  1-1 


(12) 


Equation  (12)  can  be  represented  in  matrix  notation  as  shoun  in 
Figure  2.  Fj  is  called  the  field  transfer  matrix  which  ,  by 
taking  the  limit  as  m  approaches  zero,  reduces  to  the  well- 
known  transfer  matrix  for  an  elastic,  massless  field  as  shown 
in  Figure  3. 
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Free  Body  Diagram  for  the  Feint  Mass  m^ 

Figure  4 
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Transfer  Matrix  for  a  Point  Mass 
Figure  5 
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Uith  reference  to  the  free  body  diagram  for  the  point  mass 
nrij  situated  at  point  (i)  as  shown  in  Figure  4  ,  the  state  vector 
components  on  either  side  of  rrij  are  related  as  follows: 

I  I 

M?=Mi 

N^*V>+mjWiW*  (13) 

where  w  is  the  harmonic  circular  frequency.  Thus  a  relationship 
similar  to  (12)  may  be  written 

Zf*  PjZ^  (14) 

where  Pj  is  the  point  transfer  matrix.  Equation  (14)  may  be 
written  in  matrix  notation  as  shown  in  Figure  5. 


Field 

Load 


Point 


m+n 


m+n+l 


General,  Piece-wise  Continuous  Beam 
Figure  G 

The  relation  between  the  state  vectors  Z„,^, and  Z,  at 
each  end  of  .  the  generalized  beam  shown  in  Figure  6,  having 
ni  uniform  elastic  fields  and  n  point  loads  will  be  of  the  form 


Z 


=  P,-F,-P,-F,-F,-F,- 
*U-Z, 


•Fn,-P„-Z, 


(15) 
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uhere  the  intermediate  state  vectors  will  have  been  eliminated  by 
multiplication  of  the  transfer  matrices  Pj  and  Fj  to  form  the 
overall  transfer  matrix  U  •  For  the  idealized  end  conditions 
treated  in  this  report  two  columns  of  U  can  be  ignored,  owing 
to  the  two  zero  components  ofZj;  also  two  components  of  Z^+n*!  will 
he  zero,  allowing  the  formation  of  a  2x2  determinant  D. 

A  simple  beam  may  be  treated  algebraically,  in  which  case 
the  coefficients  of  D  would  involve  functions  of  a»  .  Solutions 
for  ci)„  would  be  obtained  from  the  eigenfrequency  condition  D=0. 
Realistic  problems  demand  a  numerical  approach  .in  which  ojn  is 
surmised  and  the  value  of  D  obtained.  If  the  latter  is  plotted 
against  oi  ,  the  eigenvalues  can  be  bracketed  with  reasonable  accu¬ 
racy.  However,  using  only  eight  digits  may  give  rise  to  pre¬ 
cision  problems  at  low  frequencies  even  for  simple  beams.  A  more 
generally  applicable  numerical  approach  has  been  adapted  from 
Pestel  and  Leckie,  (ref,  1,  pp.  204-213) ,  referred  to  therein  as 
the  "modified  transfer  matrix  method" .  This  method  avoids  evalu¬ 
ating  U  and  hence  D  by  directly  obtaining  the  two  (ideally)  zero 
components  of  ,  One  of  the  two  nonzero  components  of  Z, 
is  normalized  to  unity,  the  other  being  arbitrarily  chosen.  For 
a  given  estimate  of  eigenfrequency,  a  single-valued  frequency  func¬ 
tion  R(ci))  can  be  obtained  after  two  iterations.  The  first  iter¬ 
ation  establishes  the  correction  to  be  applied  to  the  non-nor- 
malized  component  ofZ]  as  the  arithmetic  mean  of  the  two  ideally 
zero  components  of  Z,n«n*i  •  These  latter  components  are  again  ob¬ 
tained  and,  at  the  eigenfrequency  condition,  will  be  zero;  gener¬ 
ally,  however,  their  difference,  which  is  termed  R(w)  ,  will 
be  nonzero  and  can  be  considered  analogous  to  the  value  of  D  at 
a  general  <u  .  This  method  has  the  advantage  that  numerical  diffi¬ 
culties  arising  from  evaluating  D,  ,  where  the  difference  is 
several  orders  of  magnitude  smaller  than  the  minuend  and  subtra¬ 
hend,  have  been  eliminated.  This  process  was  the  principal  error 
source,  and,  using  the  modified.msthDd,  high  eigenfrequencies  of 
complex  beams  can.  be  determined  with  reasonable  accuracy  with¬ 
out  recourse  to  the  double  precision  facility. 
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Conditianing  of  Input  Data 


In  its  tirriE-sharing  form,  as  mentinned  in  the  Introduction, 
FAMSIMUB  appears  in  5  separate  parts  or  saved  items:  EMBIO  uihich 
conditions  the  input  data;  EMB15  uhich  in  conjunction  uith  EMB20 
and  EMB21  computes  the  eigenfrequencies  and  state  vectors;  and 
EMB25  which  computes  the  eigenmodes  and  stress  values.  The  fol¬ 
lowing  description  of  FAMSIMUB  is  made  with  reference  to  these 
time-sharing  component  items,  which  will  in  turn  be ’ illustrated 
bv  a  sample  problem.  This  problem  very  nearly  exhausted  the  core 
storage  capacity  of  the  Philco-Ford  system  at  the  time  of  uiritipg. 

The  beam  dealt  with  is  shown  in  Figure  7  opposite.  Seven 
parameters  PAR(I ,J) ,1=1,7  are  considered  to  describe  fully  the 
discrete  and  distributed  loading  systems  and  the  beam's  material 
and  geometrical  properties  as  follows: 


PAR(1, J) , J=l,4  are  point  weights  of  7,3,4  and  6  pounds  posi¬ 
tioned  5.1,11.9,21,8  and  32.1  inches  respect¬ 
ively  from  the  left-hand  support. 

PAR(2, J) ,3=1,2  are  uniformly-distributed  loadings  of  2  and 
3  Ib/in.,  applied  over  the  lengths  shown. 


PAR(3.J).3=1.4 
PAR(4.J). 8=1,4 
PAR(5,J).J=1.3 

PARCS, 8), 8=1, 3 
PAR(7.8).8=1,3 


are  the  values  of  material  density. 

are  the  values  of  the  elastic  modulus. 

are  the  values  of  area  moment  about  the  neu¬ 
tral  axis. 

are  the  values  of  cross-sectional  area. 

are  the  distances  from  the  neutral  axes  of 
the  various  sections  to  the  outermost  bending 
fiber. 
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Figure  9 
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The  I\I0SUP=10  different  support  conditions  are  illustrated 
in  Figure  B,  each  condition  corresponding  to  a  value  of  the  sub¬ 
scripted  variable  NSUPP.  The  state  vectors  Z  and  Zm*n*i  are 
given  in  forms  corresponding  to  the  requirements  of  the  modified 
transfer  matrix  method ,  i.e.thetuo  non-zero  components  of  Z,  ap¬ 
pear  as  unity  and  X  ,  uhich  represents  the  initial,  arbitrary  choice 
of  slope  or  bending  moment  as  the  case  may  be.  For  the  sample 
problem,  three  support  conditions  are  chosen:  IMSUPP(1)=2  (free- 
guided),  NSLIPP(2)=7  (hinged-hinged)  and  NSUPP(3)=lD  (clamped- 
clamped).  The  order  of  support  condition  can  in  general  be:  ac¬ 
counted  for  by  taking  the  origin  of  x  at  the  appropriate  end; 
thus  if  for  the  sample  problem  the  initial  support  condition 
had  been  guided-free,  the  right-hand  end  mould  have  been  taken 
as  the  origin.  It  mould  be  feasible  for'all  ten  conditions  to 
be  included  in  one  run  on  a  batch  system,  but  for  time-sharing 
purposes  the  saved  items  EMB15/EMB20/EMB21  mill  terminate  if  the 
eigenmodes  and  stress  values  are  desired,  and  EMB25  mill  have  to 
be  run  separately.  The  input  data  are  contained  in  file  EMBIDR 
as  shomn  in  Figure  9.  EMBIO  is  shomn  in  Figure  10. 

EXPLAIMATIOIM 


Line  100  mill  almays  contain  14 
items,  regardless  of  problem  com¬ 
plexity,  as  folloms: 

NSS  is  the  number  of  subspans 
into  mhich  the  beam  is  naturally, 
subdivided  by  the  point  loads 
and  the  discontinuities  intro¬ 
duced  by  the  distributed  quan¬ 
tities.  This  is  illustrated  in 
Figure  7,  mhere  the  finer,  ver¬ 
tical  lines  denote  the  inter¬ 
vals. 

IMOSUP  is  the  number  of  support 
conditions  to  be  examined. 

IMMIIMT  is  the  number  of  modes  for 
uhich  the  eigenfrequency  iB  de¬ 
sired. 


EMBIOR  Read  by 

Lines  EMBIO  Lines 

100  4-5 
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NMSHPC  <  IMMINT)  is  the  number 
of  medes  for  uhich  the  shspe 
and  stress  values  are  desired. 

MOINTP  is  a  limiter  for  the  num¬ 
ber  of  interpolations  to  be  per¬ 
formed  in  refining  the  eigen- 
frequency  value  (=10). 

G  is  the  gravitational  constant. 

BMLGTH  is  the  total  beam  length, 
to  be  checked  against  SSUM,  the 
sum  of  the  NBS  subspans. 

NPAR(I),I=1,7  are  the  numbers 
of  the  different  values  of  the 
seven  parameters  PAR(I,0)  des¬ 
cribed  above.  Thus  there  are 
1\1PAR(1)=4  different  values  of 
point  uieight,  etc. 

105-115  10-11  SS(I) ,1=1 jMBS  are  the  lengths 

of  the  subdivisions  shown  in 
Figure  7. 

120  14-16  Line  120  will  contain  a  minimum 

of  3  items  up  to  an  indeterminate 
maximum,  depending  on  the  program 
length  limitations,  as  follows::. 
NBUPP(I) ,I=1,N0BUP,  are  integer 
representations  of  the  support 
conditions  to  be  investigated 
(refer  to  Figure  8). 

MBHP( I)  ,  1=1  ,I\IMBHP  are  the  modes 
for  which  the  deflection  and 
stress  values  are  desired,  in 
addition  to  eigenfrequency.  (If 
no  eigenmodes  are  desired,  then 
MBHP(1=I\1MBHP)=0). 

125-250  23-27  The  loading,  material  and  geo¬ 

metrical  properties  as  shown  in 
Figure  7  appear  in  single-column 
format  together  with  the  corres¬ 
ponding  value  of  IPAR(I,J),  which 
represents  the  number  of  the  sub¬ 
division  at  which  the  property 
value  occurs.  Thus,  the  3  lb. 
point  load  (=PAR(1, 2))  appears  at 
subdivision  7  (=IPAR(1,2)) . 
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To  verify  that  all  input  values  have  been  read  in  from  EMBIOR 
properly,  lines  6,  12,  17,  21  and  26  of  EMBIO  create  the  file 
MBlOLdl  as  an  echo  check  (Figure  12).  EMBIO  uill  be  discussed 
ujith  reference  to  the  flow  chart  shown  in  Figure  11 ,  which  illus¬ 
trates  only  the  more  involved  steps  in  detail. 

An  error  message  (statement  153)  will  be  printed  out  and 
operations  will  be  terminated  if  SSUM,  the  total  of  SS(I)  ,I=1,IMSS 
is  not  equal  to  the  value  of  BMLGTH  as  contained  in  EMBIOR.  This 
will  afford  an  extra  check  on  the  values  of  S5(I)  as  stared  in 
core  during  execution. 

Statement  3  determines  if  a  point  l.oad  is  present  at  the 
left-hand  support,  i.e.  if  IPAR(1,1)=1.  If  not,  statements  3-6 
add  a  zero  mass  at  this  point  and  shift  the  subscript  0  of  PAR(I,0) 
and  IPARdjO)  up  by  unity.  Lines  62-66  similarly  detect  a  zero 
load  at  the  right-hand  support,  necessitating  PAR(1,I\1PAR1)=D. 

The  DO  loops  involving  lines  67-85  distribute  the  continuous 
properties  PAR(I,0) ,1=2,7  over  the  whole  beam,  giving  MSS  values 
of  each  property  (refer  to  line  69).  If ,  say ,  two  or  more  different 
values  of  uniformly  distributed  load  exist,  the  YE5  branch  of  line 
70  is  taken  in  order  to  assign  the  appropriate  value  of  3  such 
that  PAR(2,0)  is  repositioned  at  the  correct  subspan.  Thus,  for 
the  existing  problem,  PAR(2,4  )(=3  Ib/in)  ^  read  in  must  become 
PAR(2,10)  and  so  forth.  Subsequently  these  redistributed  values 
must  be  spread  out  over  the  intermediate  subspans,  e.g.  PAR(2,10) 
will  be  applied  to eubspans  11-14.  This  is  accomplished  by  lines 
77-85.  Referring  to  file  MB10LJ2  as  shown  in  Figure  13,  the  pro¬ 
perty  values  as  redistributed  can  be  seen;  also  note  that  17  (=MSS) 
values  exist  for  each  continuous  property  (1=2,7). 

.  In  order  to  establish  an  initial  guess^  to  commence  eigenfre- 
quency  computation  in  EMB15,  an  extreme  range  of  is  found  in 
line  102.  The  lower  bound, Y ,  is  based  on  minimum  stiffness  and 
maximum  distributed  loading  taken  as  continuous  over  the  entire 
span  and  using  a  guided-free  constraint.  The  upper  bound,  X  ,  is 
based  on  maximum  stiffness  and  minimum  distributed  loading,  using 
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the  clamped-clamped  constraint.  These  values  are  entered  in  file 
MBlDliJ3  as  shoun  in  Figure  14.  Their  natural  logarithms  are  desig¬ 
nated  as  FINCR,  i.B.  possible  values  for  the  initial  increment. 
The  employment  of  these  values  is  completely  discretionary. 


E'lE  1  0 


13:53  FRl.  09/05/69 


1  SPILE  EMD 1 aR,MD10Wl,MB10W2,MBI0W3 

?,  DIMENS12i\’  [.'PARC?) ,  PAR  (7,20),  IPARC  7,20)  ,PARMAX(6),PARMI  N  (6),  MSHP(20) 
+  ITEMP(20),PTEriP<20),SS(17),{1i:':T(20),NSUPP(10) 

/:  READ  (1,150)  NSS,  N3SUP,  NMINT,  MKSHP,  N0INTP,  G,BMLGTH,  (NPAR(I) ,  1=  1 , 7) 

5  150  F3RMAT(4X,5I5,2FS.3,7I3) 

6  V.'RITE  (2,100)  riSS,MOSUP,NMINT,NMSKP,N0INTP,G,BMLGTH,  (NPAR(I),I=1,7) 

7  160  F3RMAT  (291!  ECH0  CHECK  0F  INPUT  DATA/6X,4HNSS=, 13, 3X, 

8  +  6HN0SUP=  ,I3,3X,SHilMINT=,  I3,3X.6KNMSHP=  ,I3,3X,7HN0INTP=,I3/6X, 

9  +  2HG=,F3.3,3X, 7HEnL3TH=,F8.3, 13H  NPAR  VALUES, 713) 

10  READ  (1,151)  (SE(I) , 1= 1 ,NSS) 

-t  1-4  54-  -FDfmA4K4X  ,-^4r0  .-3-) - - _  _ 

12  WRITE  (2,163)  (SS(  I ) ,  1=  1 ,  NSS) 

I5  163  F3RMAT(/6X,SHSS  VALUES /( 6X, 6F 1 0. 3 ) ) 

14  READ(1,152)  (NSUPP ( I ) , 1= 1 , N0SUP) , (MI NT (I ) , 1= 1 , NMINT) , 

15  H-  (MSHPd)  ,I=1,NMSHP) 

16  152  F3RMAT(4X,  1614) 

17  WRITE  (2,164)  (f;3SUP,  (NSUPPd  ) ,  I- 1 ,  N0SUP)  )  , 

13  +  (NMINT,  (MINT(I),I= l,NMINT)), (NMSHP, (MSHP ( I ) , 1= 1 , NMSHP) ) 

19  164  FORMAT  (/6X, 12HNSUPP  VALUES/6X,*I4//6X. 1 IHMINT  VALUES 

20  +  /6X,*I4//6X, 1 IHMSHP  VALUES/6X.*I4 ) 

21  WRITE  (2, 168) 

22  163  F3RMAT(/6X, 19HPAR  AND  IPAR  VALUES) 

23  D3  200  1=  1 , 7 

24  NPR=NPAR(I) 

25  READ  (1,154)  (PAR (  I ,  J)  ,  IPARd ,  J)  ,  J=  1 ,  NPR) 

26  WRITE  (2,165)  (PAR (I , J ) , IPAR( I , J) , J= 1 , NPR ) 

27  200  CONTINUE 

23  154  F3RMAT(4X, El  1.4,13) 

29  165  F3RMA'!-(4(E1  1.4,13))) 

30  ENDFILE  2 

31  WSUM=0. 

32  NPAR1  =  NPAR(1) 

33  D3  121  1=1,  NPAR  1 

54  121  WSUM=WSUM+PAR(1 ,1) 

55  D3  122  1=2,6 

36  MPARX=NPAR(I) 

37  PAF.MAX(I)=PAR(I,  1) 

58  PARMIN(I)=PAR(I, 1) 

59  IF  (IPARd,  NPAP.X)-!)  12,122,12 

40  12  D3  122  J=2,NPARX 

41  PARMIN(I)  =  AMIHl(PARMIN(I),PARd,J)) 

42  1  22  PARMAX(I)  =  AMAX1  <PARMAXd),PARd,  J) ) 

43  SSUM  =  0i 

44  D3  1  1= 1 ,NSS 

45  1  SSUM=SSUM+SSCI) 

46  IF(ABS (BMLGTH-SSUM) - .0001)3 ,3 ,2 

47  2  PRINT  153,-SSUM 

48  155  FORMAT  (54H  INEQUALITY  0F  BMLGTH  AND  COMPUTED  BEAMLENGTH, 

49  +  ,LATTER=,F10.3) 

50  GOTO  99 

51  3  IF(IPAR(1,1)-1)4,7.4 

52  4  NPAR(l)=NPAR(l)+l 

53  NPAR  1=  NPAR  (1)  Program  vCMBlD 

» 

Figure  10 
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,  IP  CD- PAH  (1,1-1) 
t::;ip(d  =  ir/.i’c  1 ,  i- 1 ) 

PTPilPC  1  )r'' 

IT PilPC  !  ) 

D  o  I  -  1  ,  l\-  J  ill.;, 

PAKCi ,i)-PTnnP(i) 

G  IPAP  CI,I)DTEI1P(I) 

7  irCIPAPCl  ,tlPAHl)-Nr>S)8.8, 
.8  ;]PAi;ci )-iiPAnci)+i 
llPAuDllPARCl) 

•IPAnCl  .iiPARD-nCS+l 
PAHCl  ,iiPAni)  =  0. 

9  03  120  1=2,7 
|.]PAriX=IJPAR(I) 

RPAn(I)=IlS2 
IF  (NPARX-1)  11.11» 

131  IDliPARX 
1)3  10  J=2,MPARX 
IPARX=IPAR(I,K) 

PARCI  ,IPARX)  =  PAR(I,K) 

10  K=K-1 
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131 


73  10  K=K-1 

7C  11  IPARCI  ,nPARX+l)  =  NSS+l 

77-l)Tt-  A-mC - “ 

TJ  ISEP=IPAn<I  ,J-l-l)-IPAR(I,J) 
IF(ISEP-1)120, 120,135 
135  Li:i=IGEP-l 
1PARX=IPARCI,J) 

DO  120  L= 1 ,Lin 

.  «  «  T  n  A  n  /  T  I  N  Ul 

IPARX) 


70 

BP. 


DO  120  L= 1 ,Lin 
TPARY=IPAn(I , j)+L 
.  PAn(i,iPAnY)i:PAn(i, 
i  120  CONTiijiii: 

)  WRITE  (4,101)  (MPARCI) ,1= I ,7) 

»  IGl  FORdAT  (0X,11HHPAR  VALUES,  713) 

WRITE  C4,15G)  (IIPAR  1 ,  (IPARC I ,  I ) ,  1=  1 ,  NPAR 1 ) ) 
156  F3RMATC/6X, IIHIPAR  VALUES/7X, *14 ) 

WRITE  (4,157)  WSUn,GSUM 


WRITE  (4,157)  WSUn,GSUM 

157  FORllAT  (/6X,5HWGUM=,E15.8,7H  SSUM=,E15.8) 

WRITE  (3,169);  169  F0RMAT(5X,2O11PR0CESSED  PAR  VALUES) 

D 3  201  1= 1,7 
RPRrijPARd) 

201  WRITE  (5,155)  CPARd ,  J) .  J=  1 ,  NPR) 

155  FJRnAT(6X,4El5.8) 

EMDFILE  3 

EERH  3 1  =  (US  Un  /SSUIH-PARMI N  (2  )+PARf1I  N  (3  )  *PARMI  N  (6 ) )  /G 
E PR  11 02  =  ( WC  U  M  /G S  Ufl+P AR  flA  X (  2  ) +P AR 11 A X  (  5  )  * P AR M A X  (  6  )  )  /  G 
)  C3RST1=  CPARriAX(4)*PARMAX(5)/(EQRH0l*SSUM**4))**.5 
[  cailST2=CPAnt1Ill(4)*PARt1IN(5)/(EaRH02*SSUM**4))**.5 
J  WRITEC4, 167)  Xs4. 73^2*C0NST  I ,  Ys  1 . 5 71 '‘2*C0NST2 ,L0G(X)  ,L0GCY) 
J  167  r:)RnAT(/54H  EXTREME  RANGE  0F  IST/2MD  n0DE  CIRCULAR 

1  +  FREaUE:iCY/2E15.8/6HFINCR=,F7.4,llH  0R  FINCR=,F7.4) 

)  EHDFILE  4 
;  99  CONTINUE 
7  SOFT 
5  EDO 


siz: 


3300  CHARS. 


2 


NPAR(1)=NPAR(1)+1 


/  on  5  \ 
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PBIOWI 


17:42  WED.  09/03/69 


ECH0  CHECK  0F  INPUT  DATA 

NSS=  17  N0SUP=  3  NM1NT=  5  NMSHP=  5  N0INTP=  10 

G=  386.000  BMLGTHr  36.000  NPAR  VALUES  45  4  43  3  3 

SS  VALUES 

3.400  1.700  2.100  2.400  1.400  .900 

3.000  3.400  .800  2.700  2.100  2.100 

1.100  3.100  1.900  1.700  2.200  _ 

NSUPP  VALUES 
2  7  10 

MINT  VALUES 

1  2  3  4  5 

MSHP  VALUES 

1  2  3  4  5 

PAR  AND  IPAR  VALUES 

.7000E+01  3  .3000E+01  7  .4000E+01  11  .6000E+01  16 

.OOOOE+Ol  1  .2000E+01  2  .OOOOE+01  5  .3000E+01  10 

.OOOOE+01  15 

.  lOOOE+00  1  .3000E+00  4  .2000E+00  9  .5000E-01  17 

.3000E+05  1  -3500E+08  6  .lOOOE+08  12  .3300E+08  14 

.4500E-03  I  .5210E-02  8  .0125E-03  13 

.6000E-01  I  .2500E+00  8  .1500E-01  13 

.  1500E+00  1  .2500E+00  8  .5000E-01  13 


Permanent  File  MBIOLJI 
Figure  12 
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13;35 


THUR.  0S/U/S9 


?r,:cz;ssi:D  values 

.  J  . OOOOOCQOE+01  . 70n00000£+01,  . 30000000E+0 1  . 40000000E+0 1 

.SOOOQO'OOE+01  .OOOOOOOOE+01 

I.00000300E+01  .20DC0000E+01  . 20000000E+0 1  . 20000000E+0 1 

.OOOOOOOOE+01  ..OOOOOOOOE+01  .  OOOOOOOOE+0 1  OOOOOOOOE+0 1 
.OOOCOOOOE+01  .309C0000E+01  . 3 OOOOOOOE+0 1  .30000000E+0 1 

.30000000E+01  .30000000E+01  . OOOOOOOOE+0 1  . OOOOOOOOE+0 1 

.OOOOOOOOE+01 

I.lOOOOOGCE+00  . lOOOOOOOE+00  . 1 OOOOOOOE+00  .30000000E+00 
.30000C00E+00  .30000000E+00  .30000000E+00  . 30000000E+00 

.2COOOQ0OE+OO  . 20000000E+00  .20000000E+00  .20000000E+00 

.20000C00E+0C  .20000300E+00  . 2 0000 00 OE+0 0  . 20000000E+00 

.OOOOOOOOE-01 

I.30000000E+OS  .30000000E+08  . 30 OOOOOOE+OS  . 30000000E+0S 

.30000COOE+OC  .35000000E+08  . 55 OOOOOOE+08 '  .35000000E+08 

.35000000E+08  .350000C0E+08  . 35 OOOOOOE+08  . . 1 OOOOOOOE+08 

. 1000000CE+08  .33000000E+08  .33000000E+08  .33000000E+08 

.33000000E+08 

I.45000000E-03  . 45 OOOOOOE- 03  . 45 OOOOOOE-03  . 45000000E-03 

.45  000000E-03  . 45 OOOOOOE-05  .  45000000E-03  . 52 1 OOOOOE-02 

.521OOC0OE-O2  . 52 1 OOOOOE- 02  .521 OOOOOE-02  . 52 1 OOOOOE- 02 

.01250000E~03  . 01250000E-03  . 0 125 000 OE- 03  . 0 125 0000E-03 

.01250000E-03 

.60000000E-01  .60000000E-01  .  . SOOOOOOOE-01  . SOOOOOOOE-0 1 

.60000000E-01  .60000000E-01  . 60000000E- 0 1  .25000000E+00 

PAR(G,J)  .25000000E+00  . 25 OOOOOOE+00  , 25 OOOOOOE+00  .25000000E+00 

.  15000000E-01  .  15000000E-01  , 15 OOOOOOE- 0 1  . 15 OOOOOOE-O 1 

.  15000000E-01 

.  15000000E+00  . I5000000E+00  .  15000000E+00  . 15 OOOOOOE+00 

. 15000000E+00  . 15C00000E+00  . 15000000E+00  .25000000E+00 

PAR(7,J):  ,25000000E+00  . 25 OOOOOOE+00  .25 OOOOOOE+00  . 25 OOOOOOE+00 

.50000000E-01  .50000000E-01  . 5 OOOOOOOE-0 1  . 5 0000000 E- 0 1 

.50000000E-01  Permanent  File  MB lD(iJ2 

Figure  13 

l•laV3  17144  VCD.  05/03/69. 

NPAR  VALUES  6  17  17  17  17  17  17 

IPAR  VALUES 

1  3  7  11  16  18 

WSUM=  .20000000E+02  SSUM=  > .36000000E+02 

EXTREME  RAN6E  0F  1ST/2ND  MODE  CIRCULAR  FREQUENCY 
. 1941811 7E+03  .21953789E+00 

FI NCRs  5.2688  0R  FI»CR=-1.5162 

Permanent  File  MB101J3 

Figure  14  , 


>S 


(This  page  left  intentionally  blank) 
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Computation  of  Ei genf requencies  and  State  Uectors 


The  saved  items  EMB15/2D/21  accept  the  data  conditioned  by 
EMBIO  and  compute  the  eigenf requencies  and  state  vectors  for  each 
MIlMTd)  and  NSIJPP(II)  without  interruption  until  the  problem  is 
completed.  If  eigenmodes  or  further  refinement  of  particular 
eigenfrequencies  are  desired,  however,  execution  may  be  terminated. 
For  efficiency,  solutions  may  begin  at  any  value  of  guessed  fre¬ 
quency  and  any  mode  at  the  discretion  of  the  user.  Array  storage 
requirements  have  been  minimized  by  this  restart  feature;  thus 
when  3C  values  of  frequency  FREQ(I)  have  been  computed  and  used 
to  determine  the  corresponding  frequency  function  DIFCD  and  more 
are  desired,  execution  is  recommenced  by  the  user  with  the  index 
I  reset  to  1  and  employing  the  previous  frequency  as  FREQCl).  If 
the  current  MII\1T=MSHP(H\I) ,  i.e.  the  eigenmode  is  desired ,  the  re¬ 
quired  data  corresponding  to  the  eigenfrequenoy  as'  computed  after 
ten  interpolations  are  entered  into  file  MB15Ll)1  and  linked  file 
MB15liJ2/MB15U3  by  EMB20  and  the  message  "RUi\l  EMB25  FDR  MODAL  V/ALUES" 
is  printed. 

The  flow  chart  of  EMB15  is  shown  in  Figure  16  and  the  program 
itself  appears  as  Figure  15.  The  input  is  largely  from  permanent 
files  EMBIDR,  MBlDklB  and  MB1DIJ3.  Terminal-supplied  input  (lines 
4B-5D)  is  confined  to  values  which  implement  the  restart  feature. 


Files 

Read  bv 

EMB15  Lines 

EXPLANATIDiM 

EMBIDR 

lD-22 

Refer  to  "Donditioning  of 
Input  Data". 

MBlDliJ2 

3D-38 

Refer  to  Figure  13. 

MBlD[il3 

24-26 

This  input  furnishes  the. 
number  of  point  loads  to  be 
considered,  including  vir¬ 

tual  loads,  if  any.  ([\1PAR 
(I), 1=2, 7  are  all  =IMSS) 

27 


MB1DLJ3  40-42  The  subspans  corresponding 

to  each  point  load,  IPAR(1,0) 
are  needed  for  the  computation 
of  point  transfer  matrices. 

Initially  and  uhen  reentry  into  the  program  is  attempted  a 
message  is  printed  requesting  seven  parameters  from  the  terminal: 

1)  II_  is  the  index  of  the  support  condition  NSUPP  and  can  be 
changed  at  will  to  alloui,  for  example,  an  investigation  for  eigen¬ 
values  uithin  a  given  frequency  range  for  various  end  constraints. 

2)  &  Ij^  are  the  subscripts  of  MIIMT  and  MSHP  respectively  and  can 
be  set  independently  to  alloiu  the  determination  of  an  eigenfre- 
quency  of  interest  for  other  than  MINT(l). 

3)  IMMDDE  denotes  the  mode  currently  being  investigated. 

4)  FREQ(l)  is  the  initial  frequency  to  be  used  as  a  starting  point 
in  the  spectrum  under  examination. 

5)  FIIMCR  is  the  natural  logarithmic  increment  of  frequency  to  be 
added  to  FREQ(l)  and  may  be  either  arbitrary  or  draun  from  the 
extreme  range  of  values  as  computed  by  EMBID.  Negative  values 
must  be  avoided. 

6)  IND(I )  indicates  the  polarity  of  the  frequency  function  DIF(I) 
(=R(  w  ))  ;  i.B.  IND(I)=1  if  DIF(I)>0,  otheruise  IND(I)=2.  IND(I) 
is  needed  to  initialize  the  procedure  in  EMB20  for  detecting  po¬ 
larity  reversal,  i.e.  if  (IND(I)-IND(I-l)).  The  user-  uill  initi¬ 
ally  not  knou  the  sign  of  DIF(l)  and  thus  a  restart  may  be  indica¬ 
ted  if  the  wrong  choice  is  made. 

In  line  42  an  arithmetic  statement  function  is  defined  which 
represents  the  procedure  for  interpolating  between  the  previous 
frequency,  FREQ(I-l)  and  FREQ(ICROSS) ,  the  last  frequency  corres-r 
ponding  to  a  DIF  value  of  different  sign. 

In  lines'  52-54  the  frequency-independent  B(J),J=1,NSS  are 
computed  for  use  later  .in  establishing  field  transfer  matrices. 
Lines  58-68  set  FRQN(1)=D  for  the  first  four  support  conditions 
and  increment  IM  and  IN  to  2  (if  MINT(1)=1  through  oversight).  The 
SS(J.X  values  have  not  been  included  in  B(J)  as  an  expedient  for 
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saved  item  EMB25,  uhich  reads  the  CDEFM(J)  and  further  subdivides 
the  subspans;  thus  it  is  desirable  that  T  (=COEFM(d)=B( J)*FREC1(I )**)^) 
does  not  involve  SS(J)  at  this  stage.  If,  at  line  56,  I\1M0DE>2,  the 
immediatelv  preceding  eigenfrequency  uill  be  supplied  at  the  terminal. 
This  is  utilized  by  EMB2D  to  compute  a  FINCR  for  the  ensuing  mode, 
i.e.  for  I\1MDDE+1. 

In  line  72,  FREQ(2)  is  determined  exponentially;  arithmetic  in¬ 
crementing  mould  be  less  efficient  in  the  preliminary’  detection  of 
eigenvalues.  Care  should  be  taken  not  to  supply  FRC}N(1)=D  to  avoid 
an  execution  error  in  connection  mith  the  intrinsic  LOG  function: 

Statement  16  is  a  computed  GO  TO  mhich  initializes-  the  first 
state  column  vector  STAMAT  and  correction  column  vector  CORMAT  accor¬ 
ding  to  the  support  condition.  Thus  for  the  guided-free  constraint 
(I\!SUPP(  1)=2)  the  initial  vectors  are: 


STAMAT^  =  |l  1  0  □! 
CORMATj^  =  {□  1  □  D| 


(16) 


Here  the  free-end  deflection  is  normalized  to  one  inch  and  the  slope 
chosen  as  unity.  Lines  84-128  establish  the  transfer  matrices,  in¬ 
dices  H  and  8  being  incremented  at  each  point  load  and  elastic  field 
respectively.  Figure  24  illustrates  horn  the  K  and  3  values  are  assigned 
for  this  beam.  Llhenever  3=IPAR(l,li)  a  point  transfer  matrix  is  com¬ 
puted  in  lines  88-90.  If,  subsequently,  IMSB  <  IPAR(I,K) ,  (line  92) 
then  matrix  computation  is  completed  for  FREQd)  and  control  passes 
to  EMB21  to  obtain  the  state  vectors  and  correction  factor  by  chain 
multiplication;  othermise  computation  of  the  Jth  field  matrix  is  in¬ 
dicated — .according  to  Figure  3  if  the  field  is  effectively  massless 
(i.e.  if  PAR(3,3)  < .0001)  or  Figure  2  othermise.  The  argument , of 
the  natural  trigonometric,  and  hyperbolic  functions  assumes  its  larg- 
est  value  at  SS<13),  i.e.  ARG  =  .097724  a)  and  the  largest  per- 
missible  argument  =  }176|  ,  uhich  corresponds  to  =  517  kHz  approx¬ 
imately.  Argument  magnitude  difficulties  may  be  present  uith  beams 
'having  longer  and  more  compliant  elastic  fields.  In  such  cases  a 
virtual  load  could  be  inserted  at  midspan,  thus  halving  the  argument. 
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For  example,  the  natural  frequencies  in  pure  bending  of  a  1"  square 
uniform  aluminum  beam  25"  long  and  clamped  at  each  end,  as  reported 
in  Reference  2,  uere  corroborated  to  nearly  27  DOD  rad/sec  (until 
the  argument  exceeded  |17S1  ). 

EMB21  appears  as  Figure  17  and  is  charted  in  Figure  18.  Its 
function,  as  an  item  which  is  compiled  with  EMB2D  and  EMB25,  is  to 
obtain  the  intermediate  and  final  state  vectors  (STAMATC J,K,1) , 
J=2,NSS+NPAR(1)+1) ,H=1,4  and  the  elements  of  the  correction  vector 
CORMAT(J],h,l).  At  each  FREQCI)  as  mentioned  previously  two  itera¬ 
tions  are  needed  to  obtain  DIF(I);  ITER=D  causes  both  vectors  to 
be  computed  in  lines  154-16D.  If  ITER=1  only  STAMAT  is  recomputed. 
By  the  computed  GO  TO  in  line  162  the  components  CORl  and  C0R2 
of  the  correction  term  COR  are  obtained  in  accordance  with  the  sup¬ 
port  condition.  Thus  for  this  beam  slope  and  shear  force  at  the 
right-hand  support  are  ideally  =0,  causing  M=2  and  l\)=4.  CORl  and 
C0R2  are  obtained  from  the  relationship 

STAMAT(LIM+1 , 2 , 1)+C0R1*C0RMAT(LIM+1 , 2 , 1)=0 

~'^m+n+l  (17) 

STAMAT(LIM+1 , 4 , 1)+C0R2*C0RMAT(LIM+1 , 4 , 1)=0 

~ '^m+n+l  (IB) 

The  correction  term  COR  is  taken  as  the  arithmetic  mean  of  CORl 
and  C0R2  and  is  applied  in  statement  105  to  the  left-hand  slope, 
the  other  elements  of  remaining  the  same.  ITER  is  now  made  =1 
and  statement  300  is  reentered.  After  the  recomputation  of  STAMAT, 
exit  is  made  to  statement  119  in  either  EMB20  or  EMB25  as  appro¬ 
priate'.  In  EMB20  (refer  to  Figures  19  and  20)  statement  11^  com¬ 
putes  the  frequency  function  as  the  difference  between,  the  compo¬ 
nents  of  (which  should  be  zero  at  the  eigenf requency) .  Ac¬ 

cording  to  {i }  this  difference  is  to  be  takerf  as  C0R1-C0R2,  but 
this  involves  magnitudes  that  become  progressively  smaller  as  NMODE 
increases.  This  results  from  the  growth  of  the  CORMAT  elements 
with  IMMQDE,  causing  the  graph  of  DIF  vs.  FREQ  to  flatten  out.  In 
some-  cases  this  may  hinder  eigenfrequency  determination  if  suffi¬ 
ciently  severe. 
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statement  80  causes  central  to  branch  to  84  if  interpalation 
is  not  being  performed,  i.e.  if  the  frequency  axis  has  not  yet  been 
crossed  for  the  current  NM0DE=MI(\1T(IM) .  The  test  for  polarity  re¬ 
versal  of  DIF(I)  is  statement  84 — the  YES  branch  being  taken  to  ob¬ 
tain  additional  DIF  values  until  a  reversal  does  occur,  whereupon 
the  current  I  is  retained  in  MAX  to  preserve  II\1D(I)  and  also  allow 
frequency  incrementing  to  commence  with  FREQ(MAX)  after  all  inter¬ 
polations  are  completed.  The  first  approximation  to  FRQNClMMODE)  is 
a  linear  interpolation  between  FREQ(I)  and  FREQ(I-l)  and  is  printed 
by  lines  21  and  23.  After  this  first  axis  crossing,  if  l\IMODE=Mir\lT(IM) 
the  YES  branch  is  taken  to  97  and  ten  (=I\10I[\1TP)  further  interpola¬ 
tions  are  performed  to  refine  FRQl\i.  I-l  is  retained  in  ICROSS  to 
preserve  the  FREQ  value  previous  to  the  axis  crossing.  The  first 
approximation  to  FRQN(IMMDDE)  is  put  into  FREQ(I+1)  and  control  is 
returned  to  statement  16  in  EMB15  to  obtain  the  corresponding  DIF. 
Interpolation  proceeds  until  INTERP>N0INTP  ,  at  which  point  the  final 
value  of  FRQN  is  equated  to  the  current  FREQ.  If  a  sign  change  in 
DIF  occurs  during  interpalation  and  INTERP  <  IMOINTP ,  then  lines  37 
and  39  assign  the  current  I-l  to  ICROSS  and  control  is  returned  to 
EMB15. 

If  the  current  eigenf requency  was  the  last  to  be  determined, 
i.e.  NM0DE=MINT(NMII\1T)  ,  then  statement  133  determines  if  the  eigen- 
mode  is  also  required,  i.e.  if  l\IMODE=MBHP(NMBHP) .  The  NO  branch 
from  133  leads  to  the  question  of  whether  or  not  all  support  condi¬ 
tions  have  been  investigated,  i.e.  if  NOSUP  >  II .  Here  the  NO  branch 
indicates  problem  completion ,  a  notice  to  this  effect  being  printed 
by  statement  134;  otherwise  139  causes  a  request  to  be  printed  for 
the  new  NSUPP,  and  control  is  returned  to  input  statement  12  of 
EMB15.  The  YES  branch  from  133  causes  all  data  needed  by  EMB25 
to  be  entered  into  permanent  files  MBISUI  and  MB15U2/MB151J3  by  lines 
79-99.  This  involves  the  state  vectors  STAMAT  and  all  transfer  ma¬ 
trices  BMAT ;  format  statements  75,  76,  78  and  79  arrange  the  ele¬ 
ments  to  facilitate  visualization  of  the  matrix  multiplication  (refer 
to  Figure  21  which  shows  the  first  five  state  vectors  for  the  fourth 
eigenTrequency  of  156.7256  rad/sec  (=24.944  Hz)).  MBlSliJl  is  shown 
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in  Figure  22,  FRQN(NMODE)  is  added  for  completeness  and  FINCH  will 
be  required  when  EMB15  is  rerun. 

Lines  59-77  establish  the  new  increment  for  the  succeeding  NMODE 
on  the  basis  of  support  condition  and  the  value  of  NMODE.  Thus 
for  NM0DE=1  or  2  and  NATSP=1,2,3  or  4  there  exists  no  FRQN(NMODE-l) 
as  required  by  statement  72,  causing  FINCR=,1,  Otherwise  statement 
72  assigns  a  value  to  FINCH  based  on  the  log  ratio  of  the  current  to 
the  previous  FRQN  if  NM0DE=2;  if  NM0DE=1  for  NATSP=5-10  then  FINCR=.l. 
After  FINCR  has  been  determined  by  statement  70  or  72  it  is  recorded 
in  MB151J1  if  the  eigenmode  is  desired,  otherwise  line  73  increments 
I ,  IM  and  NMODE.  The  IND(I-l)  and  FREQ(I )  values  are  based  on  IND(MAX) 
and  FREQ(MAX). 

Figure  23  shows  the  results  for  NM0DE=NATSP=2.  The  sensitivity 
of  DIF  to  FREQ  is  quite  apparent  at  1=19  and  20,  where  9th  and  10th 
place  differences  in  FREQ  cause  an  80,2%  decrease  in  DIF. 
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(This  page  left  intentlonallv  blank) 
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1^;?4 


WED,  09/05/69 


"FILE  E^iB10R,!1B10W2,iv)B10W3,MB15Wl  ,MB15W2/MB15W5 
^  DI«1E^!SIPIN'  Bf'1AT(23,4,4),C0EFn(17),FREQ(5O)  ,DIF(30)  ,IND(30)  , 

FRQrK50) ,niNT(30) ,MSHP(30) , NSUPP(IO) ,MPAR(7) ,BS( 17)  , 

P  +  ?AR(7,2O),IPAR(1,2O),B(17),STAMAT(24,4,1),C0RMAT(24,4,1) 

10  READ  (  1,290)  MBS ,  N0SUP,  Nf1 1  NT,  NMBHP,  M0 1  MTP,  G 
12  290  FOR'VIAT  (4X,  5  15  ,  F9.3) 

14  READ  (  1,292)  (SS ( I ) , I = 1 , MSS ) 

16  292  FORMAT  (4X,6F10.3) 

15  READ  (1,291)  (NSIIPP(I),Irl,N0SUP),(^^INT(I),I=l,NMINT), 

20  +  (MSHP(I) , 1= 1 , MMSHP) 

22  291  FORflAT  (4X,  1614) 

24  READ(3,293)  ( NPAR ( I ) , I r 1 , 7) 

26  295  FORMAT  (17X,7I3) 

28  READ(2) 

30  D0  294  1=1,7 
32  MPR:M?AR(I) 

34  NPAR1=NPAR(1) 

36  294  READ  (2,295)  (PAR( I , J) , J= 1 , MPR) 

38  295  FORMAT  <6X,4E15.8) 

40  READ(3,297)  ( NPAR I , ( IPAR( 1 , I ) , I = 1 , NPAR 1 ) ) 

42  297  F0RMAT(//7X,*I4) ;  R( I ) =- ABS(DIF( I- 1 ) ) / ( ABS(DIF( I- 1 ) ) 

44  +  +ABS(DIF(ICR0SS)))*(FREQ(I“1)-FREQ(ICR0SS))+FREQ(I-1) 

46  LIM=MSS+NPAR(1);  12  PRINT, "INPUT  VALUES  OF  1 1 , IM, IN, NM0DE, FREQ( 1 ) 
48  +FINCR  AND  IND(I)"?  I NPUT, 1 1 , IM, I N , NM0DE , FREQ( 1) , FI  NCR , I ND( 1 ) 

50  NATSPrNSUPPCII);  ITER= INTERPrO;  1=1 

52  D0  11  J=1,NSS;  11  B<J)=(<PAR(2,J)+PAR(3,J)*PAR(6,J))/ 

54  +  (G*PAR<4,J)*PAR(5,J)))**.25 

56  IF(NM0DE-1)  201,201,202;  202  PRINT  203, "FRQN( " , NM0DE- 1 , ") =" 

58  205  F0RMAT<I2);  INPUT,  FRQN( NM0DE- 1 ) ;  Q0  TO  15 
60201  G0  T0<13,13,13,13,15,15,15,15,15,15),  NATSP 
62  15  FRQN(1)=0. 

64  NM0DE=2 

66  IF(MINT( 1 )-l )  15,200, 15 
68  200  IM=IN=2 
70  15  1=1+1 

72  FREg(I)=EXP(L0G<FREQ(I-l))+FINCR) 

74  16  GO  TO  (25,25,25,26,26,28,29,28,28,28),  NATSP 
76  25  STAMATd,  1,  l)rSTAMAT<l,2,l)  =  l.;  G0  TO  27 
78  26  STAMATd, l',l)=STAMAT(l, 5, 1)  =  1.;  G0  TO  52 
80  28  STAMAT(l,5,.l)  =  STAMATd,4,  l)  =  i.;  GO  TO  32 

82  29  STAMATd,2,l)=STAMATd,4,l)  =  l.;  27  CORMATd  ,2,  1 )  =  1  , ;  GO  TO  39 
84  32  CORMATd, 5,  Drl.;  39  K=l;  J  =  0;  17  J=J+1 
86  .IF(J-IPAR( 1,K))  23,20,23 

38  20  DO  22  L=l,4;  22  BMAT( J+K- 1 , L,L) = 1 . ;  BMAT(J+K-1 ,4, 1 )=PARd .K) 

90  +  ♦FREQ(I)**2/G 

92  IF(NSS-IPARd,K))  24,19,19;  19  K  =  K+1 

94  23  N.=J+K-1;  EI  =  PAR(4,J)*PAR(5,J);  IF(PAR(3, J)-.0001 )  5,5,6 
96  5  C0EFM(J)  =  O,;  D0  2  L=l,4;  D0  1  M:l,4;  1  Bi1AT(N,L,M)  =  0, 

98  2  BMAT(N,L,L)=1 .;  BMAT(N, 1 ,2)=BMAT(N,3,4)=SS(J) 

100  BMAT(N, 1 ,3)=BMAT(N,2,4)=SS(J)**2/<2.+EI);  BMATC N, 2 . 3) =SS( J) /El 
102  BMATCN, 1,4)=SS(J)**3/<6.*EI);  60  T0  17 

Saved  Item  EMB15 
Figure  15 
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104  s  T  =  CP!FF!1(.J)=B(J)*FRFQ(I)**,5;  ARG=T*SS(J) 

106  Ur  (FXP(ARG)+FXP(-ARG))/?..;  V:CP)S(ARG);  W=U-V 
lOf?  Xr  (EXP(ARG)-FXP(-ARG))/2.;  Y  =  SIN(ARG);  Z:X-Y 
110  3!^AT(M,  1  ,  l)rRMAT(M,2,?.)rB^1AT(N,3,3)=BMAT(N,4,4)r(U+V)/2. 
112  B^IATCN,  1 ,2)rBMAT(N,3,4)r(X+Y)/(2.*T) 

1 1  4  3X AT (  N ,  1 , 3 )  rBriAT (  N ,  2 , 4 )  r  W/(2 . I*T**2) 

116  BflATCM,  1  ,4)rZ/(2.*EI*T**3) 

UP  3MAT(M,2,  l)rpNlAT(N,4,3)rT>kZ/2. 

120  3MAT(\',2,3)r(X+Y)/(2.*FI*T) 

122  B:^AT(m,3,  n=RXAT(N,4,2)rEI*T**2*Z/2. 

124  3MAT(N,3,2)rEI*T*Z/2. 

126  BMAT(M,4, l)rEI*T**3*(X+Y)/2. 

122  G?  T0  17 
130  24  ^UPE  EnB2C 


LF.MGH 

AB0UT  3000  CHARS. 
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FIquj  Chart  of  EMB15 
Figure  16 


\  READ 
\  EMBICR 
\  MB10LJ2 
\mB101i13  / 


IIMPUT 

\  II  ,IM,I[M,i\IMCDE,  / 
FREQ(1),FIMCR  / 
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5_ _ 

COMPUTE  TRAIMSFER  MATRIX  ELEMENTS 
FOR  MASSLESS  ELASTIC  FIELD 
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F!1821  4:38  SAT.  09/06/69. 

150  300  m  310  J=1,LIM;  D0  310  M=l,4;  STAMATC J+1 , M, 1 )=0. 

152  IF(ITER)  310,180,310;  180  C0RMAT<J+1 ,M, 1 )=0.;  ilO  CONTINUE 
154  D0  340  J=1,LIM;  D0  340  M=l,4;  D0  340  N=l,4 
156  STAMAT(J+1 ,M, 1) =ST AMAT( J+ 1 , M, 1 )+ST AMAT ( J, N, 1 )*BMAT( J, M, N) 
158  IF(ITER)  340,350,340;  350  C0RMAT< J+ 1 , M, 1 )rC0RMAT( J+ 1 , M, 1 )+ 
160  C0RMAT(  J,  N,  1)*BMAT<J,M,  N);  340  CONTINUE  . 

162  G0  T0  (113,114.115.114,115,113,115,115,114,116),  NATSP 
164  113  M=3;  N=4;  60  T0  118 

166  114  M=2;  N=4;  G0  T0  1 18 

168  115  M=l;  N=3;  G0  T0  118 

170  116  M=l;  N=2;  118  IFCITER)  119,121,119;  121  C0R1= 

172  +  -STAMATCLIMfl.M,  l)/C0RMAT(LI«f  1,M,  1) 

174  C0R2=-STAMAT(LII1fl,N,l)/C0RMAT(LIM+l,N,l) 

176  C0R  =  <C0Rl+C0R2)/2. 

178  G0  T0  (105, 105, 105, 106, 106, 106, 105, 106, 106, 106), NATSP 

180  105  STAMATd  ,2, 1  )  =  STAMAT(1 ,2, 1  )+C0R;  G0  T0  330 

182  106  STAMATd  ,3, 1  )=STAMAT(1 ,3, 1)+C0R;  330  ITER=1;  G0  T0  300 


Saved  Item  EMB21 
Figure  17 
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Flow  Chart  of  EMB21 


39 


EMB20 


,18502  M0N.  08/18/69 


I  $USE  EMB21 

3  119  DIF<I)=STAMAT<LIM+1,M,1)-STAMAT(LIM+1,N,1)J  ITER=0 
5  PRINT  401,  "difference  VALUEC", DIF(I) FREQ=",FREQ(I) 

7  401  F0RMAT<I2,2E15.8) 

9  IFCDIFd))  81.81,83 

II  81  IND(I)=2|  Q0  T0  80;  83  IND<I)=1 
13  80  IF(INTERP)  125,84,125 

15  84  IF(IND(I-i)-IND(I))  96,15,96 

17  96  MAX=I;  FRQN<NM0DE)=ABS<DIF(I-!))/(ABS(DIF(I))+ABS(DIF(I-l))) 

19  +  *<FREQ(I)-FREQ(I-1))+FREQ<I-1) 

21  PRINT  4O2.*"R0UGH  FRQN<"-NM0DE,")=".FRQN(NM0DE) , "=  ", 

23  +  FRQN(NM0DE)/6.2832,"HZ  402  F0RMAT( 12. 2E 15 ,8) 

25  IF(NM0DE-MINT(IM))  135.97,135 

27  97  INTERP=1;  ICR0SS=I-1;  PRINT  403,  ICR0SS;  1=1+1 
29  FREQ<I)=FRQN(NM0DE);  G0  T0  16 
31  403  F0RMAT(7HICR0SS=,I4/) 

33  125  IFCN0INTP-INTERP)  131,131,126 
35  126  IF(IND<I)-IND<I-1))  130,137,130 

37  130  ICR0SS=I-1;  PRINT  4O3.ICR0SS;  157  lrI+1;  FREQ(I)=R(I) 

39  INTERP=INTERP+1;  G0  T0  16 
43  131  FRQN(Nn0DE)=FREQ(I) 

45  PRINT  405,*"FINAL  FRQN<"-NM0DE,")=",FRQN(NM0DE>,"=  " , FRQN( NM0DE) 

47  +  /6.2832,"HZ","  INTERPr  .INTERP;  405  F0RMAT<I2.2E15,8, 13//) 

49  INTERPsO 

51  IF<MINT(NMINT)-NM0DE)  133,133,135;  133  IF<NM0DE-MSHP(IN)) 132,74, 132 

53  132  IF(N0SUP-II)  134,134,139;  139  PRINT, "BEGIN  NEW  SUPP0RT 

55  +  C0NDITI0N";  G0  T0  12 

57  134  PRINT, "PR0BLEM  COMPLETED";  60  T0  304 

59  135  IF(NM0DE-2)  136,136,72 

61  136  Q0  T0  (70,70,70,70,71,71,71,71.71,71),  NATSP 
63  70  FINCR=.l 

65  IF(NM0DE-MSHP<IN))  73,74,73 
67  71  G0  T0  (7O,72),NM0DE 

69  72  FINCR=.1*L0Q<FRQN<NM0DE)/FRQN(NM0DE-1)) 

71  IF(NM0DE-MSHP<IN))  73,74,73 

73  73  1=1+1;  IM=IM+1;  IND(I-1)=IND<MAX);  NM0DE=NM0DE+1 
75  FREQ(I)=EXP(L0Q<FREQ(MAX))+FINCR) 

77  40  T4  16 

79  74  WRITE  (5,75);  75  F0RMAT(59X,9HSTAMAT( 1 ) ) 

81  WRITE  (5,76)  (STAMAT( 1 ,M, 1 ) ,M= I ,4);  76  FORMAT  (57X,E14.8) 

83  D0  77  J=1,LIM;  WRITE  (5,78)  J,J+1 

85  78  F0RMAT(/5HBMAT(,I2,1H),51X,7HSTAMAT(,I2,1H)) 

87  D0  77  K= 1,4;  WRITE  (5,79)  (BMAT(J,K,L) ,L= 1 ,4) ,STAMAT(J+1 ,K, 1 ) 

89  79  F0RMAT(4E14,8,X,E14.«);  77  CONTINUE;  ENDFILE  6 
91  WRITE(4,69)  NM0DE,MINT(NMINT) , II, NATSP, FRQN(NM0DE) , 

93  +  (C0RMAT(1,M,1).M=1,4),(J,C0EFM(J),J=1,NSS) 

95  69  F0RMAT(6HNM0DE=,I2,14H  MINT(NHINT)=. I2,5H  II=,I2. 

97  +  8H  -NATSP=,I2,7H  FRQN=,E15.8/9HC0RMAT(1) , 4F4. 1 
99  +  /(6HG0EFM(,I2.3H)s  ,E15.8));  ENDFILE  4 
101  PRINT, "RUN  EMB25  F0R  H0DAL  VALUES" 

103  304  CONTINUE 
105  $0PT  SIZE 

107  END  Saved  Item  EMB20 

Figure  19 
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FIdui  Chart  of  EMB2D 
Figure  20 
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FIiMCR  =  .l*ln 


FRQMdMMDDE) 


FRQN(MMDDE-l) 


(  NMDDE^MSHP(II^J))^ 


no 
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UJRITE 
STAMAT(J+l,Ki,l)' 
&  BMAT(J,H,L)  IN 
MB15uJ2/MB15Ld3 


‘  LiJRITE 

NMDDE,MII\1T(NMINT)  , 

II  ,NATSP,FRQI\1(NMDDE)  , 
]ORMATCl,M,l),CDEFM(B), 
IN  MB15ljJl 


PRINT ' 
"RUN  EMB23  FDR  MODAL  UALUES" 


END 
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!»B15W1 

NM0DE=  4 
C0RMAT<1) 
C0EFM(  1) 
C0EFM(  2) 
C0EFM(  3) 
C0EFMC  4) 
C0EFM(  5) 
C0EFM<  6) 
C0EFM<  7) 
C0EFM(  8) 
C0EFM(  9) 
C0EFM(1O) 
C0EFM(11) 
C0EFM( 12) 
C0EFM<13) 
C0EFM(14) 
C0EFM(15) 
C0EFM(lS) 
C0EFM(17) 


10:45  TUE.  09/23/69. 

MINT(NMINT)=  5  11=  1  NATSP=  2  FRQN=  . 1567256 1 S+03 

0.0  1.0  0.0  0.0 
.72925 182E-01 
.31 18332 7E+00 
.31183327E+00 
.3 122985 7E+00 
.95974937E-01 
.92346649E-01 
.92346649E-01 
.71525685E-01 
.64630752E-01 
. 18062243E+00 
. 18062243E+00 
.24705218E+00 
.1 1119480E+01 
.82500390E+00 
.14667209E+00 
. 14667209E+00 
. 10371283E+00 


Permanent  File  MB15U1 
Figure  22 
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RUM 


WED.  09/05/69. 


rj  E'1320 
I.M  EMB21 
IM  F'^B20 
SIZE  AT  LTME 


105,  3009,  1^10, 


INPUT  VALUE 
?  ?  1  9  9  '> 
FRQM(  l)r 
?  ?0. 

DIFFERENCE 

DIFFERENCE 

DIFFERENCE 

DIFFERENCE 

DIFFERENCE 

DIFFERENCE 

difference 

DIFFERENCE 

DIFFERENCE 


S  CF  II,If1, 


IN,NM0DE,FREQ(1),FINCR  AND  IND(l) 


VALIJEC  2): 
VALIJE(  3)r 
VALUE(  4)r 
VALUE (  5): 
VALIIE(  6)  = 
VALUEC  7)= 
VALUEC  8)= 
VALIJEC  9)  = 
VALUEC 10)r  - 


.2175271RE+01  FREQ: 
.31549069E+01  FREQr 
.43915608E+01  FREQ: 
.  57-98224  lE+01  FREQ: 
.7081 1595E+01  FREQ: 
.75650462E+01  FREQ: 
.59752824E+01  FREQ: 
.35757615E+00  FREQ: 
.1 1560836E+02  FREQ: 


.6lb70138E+01 
. 74591 255E+01 
.91 105939E+01 
. 1112770 5E+Q2 
. 13591 409E+02 
. 16600584E+02 
.20275999E+02 
.24765162E+02 
.30248236E+02 


ROUGH  FRQNC  2):  .24929665E+02=  . 5967670 1 E+0 1  HZ 


ICR0SS:  9 

DIFFERENCE  VALUEClDr  .75885236E-0 1  FREQ:  .24929665E+02 
ICR0SS:  10 

DIFFERENCE  VALUEC 12):  . I 5826825E-0 1  FREQ:  .24964348E+02 

DIFFERENCE  VALUEC13):  .32955923E-02  FREQ:  .2497 1 572E+02 
DIFFERENCE  VALUEC 14):,  .66740 1 58E-05  FREQ:  .24973075E+02 

DIFFERENCE  VALUEC 15):  . 1 77 1 0266E-05  FREQ:  .24973380E+02 

DIFFERENCE  VALUEC 16):  - .22753404E-04  FREQ:  .24973461E+02 
ICR0SS:  15 

DIFFERENCE  VALUEC 17):  .47282049E-04  FREQ:  .24973452E+02 

ICR0SS:  16 

DIFFERENCE  VALUEC18):  . i 4765 1 4 1 E>04  FREQ:  .24975458E+02 
DIFFERENCE  VALUEC 19):  .30257404E-05  FREQ:  . 24975459E+02 

DIFFERENCE  VALUEC20):  .59870047E-06  FREQ:-  ,24973459E+02 

FINAL  FRQNC  2):  .24973459E+02:  . 39746402E+0 1  HZ  INTER?:  10 

RUN  EMB25  FSL  Ti^DAL  VALUES 

Sample' Output  for  Secand  Mode 
Figure  23 
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(This  page  left  intentionally  blank) 
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Computation  of  Eigenmodes  and  Stress  V/alues 

Lines  1-57  represent  the  input  section  of  EMB25  (refer  to 
Figures  25  and  26),  uhich  is  based  on  the  initial  data  set  EMBIDR, 
the  output  files  MBlObJE  and  MB1DIjJ3  of  EMBIO  and  the  output  files  of 
EMB15/20.  The  program  computes  (UAL (I  ,L)  ,1=1,6)  ,L=1  ,LSUM+i\ISS+2)  re¬ 
presenting  the  four  components  of  the  state  vectors  plus  the  maximum 
bending  and  average  shear  stresses  at  each  of  the  53  values  of  dis¬ 
tance  DIST(L)  from  the  left-hand  support.  The  beam  as  partitioned 
by  the  point  loads  and  property  jump  discontinuities  is  shoujn  in 
Figure  24,  mhich  also  shams  the  53  stations  corresponding  to  the  L 
values.  At  each  discontinuity,  L  is  incremented  to  account  for 
shear  force  change  due  to  the  point  load  as  mell  as  for  stress 
changes  due  to  the  possibility  of  variations  In  cross-section  pro¬ 
perties  PAR(5-7,B).  The  latter  occur  at  L=22,23  and  L=37,38;  the 
remaining  1C  discontinuities  do  not  affect  stress  values  as  they 
correspond  to  changes  in  distributed  load,  material  density  or  elas¬ 
tic  modulus  PAR(2-4,J).  UAL(I,L)  are  first  computed  on  the  basis 
of  existing  state  vectors,  i.e,  the  STAMAT  values  as  read  in.  The 
largest  deflection  (in  absolute  value)  is  then  found  (UALMAX)  and 
compared  uith  unity  to  form  RATIO.  The  appropriate  normalized  com¬ 
ponent  of  Z^,  either  liJ  or  U,  is  then  multiplied  by  RATIO  and  a  neui 
series  of  state  vectors  obtained  through  the  $USE  EMB21  monitor 
line.  These  renormalized  values  are  returned  to  statement  119  and 
a  neu  sequence  of  UAL(I,L)  obtained  uhich  mill  contain  a  maximum 
displacement  of  +  1  inch. 

Lines  59-69  establish  the  number  and  length  of  the  subspan 
segments,  XI\IOBEG(J)  and  SUBSEG(B)  respectively;  a  maximum  of  20 
segments  is  assured  (line  50),  but  this  may  be  varied  by  the  user 
at  mill.  For  the  given  beam  no  subspan  > SSUM/20 ,  hence  XIMOSEG(J) 
=2-  for  all  B.  Line  71  establishes  the  limit  of  L,  LIM2,  as  the 
sum  of  the  number  of  subspans  NSS,  the  total  number  of  segments  LSUM, 
plus  one  station  for  each  end  of  the  beam.  Lines  73-75  initialize 
UAL(I,L)  using  the.  state  vector  as  read  in  or  returned  by  EMB21. 
In  line  61  J=IPAR(1,K)  means  that  UAL(I,L)  are  to  be  computed  for 
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the  right  side  of  a  point  load,  and  if  I\ISS  <  IPAR(1,K)  lines  85-87 
compute  the  UAL(I,LIM2)  and  control  passes  to  statement  103.  If 
0i^IPAR(l,h)  lines  93-97  compute  the  UAL(I,L)  for  the  left  side  of 
a  general  subspan.  Lines  99-147  are  involved  in  the  computation  of 
UALdjL)  for  all  the  remaining  stations  mithin  a  subspan.  SPAN  is 
a  summer  of  the  segment  lengths  for  each  subspan  and  has  a  maximum 
of  SS(J).  The  value  of  L  at  the  beginning  of  a  subspan  is  retained 
in  LSAV/E,  to  which  the  current  SPAN  is  added  to  obatia  the  corres¬ 
ponding  DIST(L)  (line  103).  As  long  as  M<  LIMl  control  passes  to 
statement  44  which  ccmsiders  the  possibility  of  a  massless  field. 
The  NO  branch  causes  V/AL(I,L)  to  be  computed  in  accordance,  with 
Figure  3 — the  YES  branch  in  accordance  with  Figure  2;  control  then 
passes  to  42.  When  M=LIM1  the  V/AL(I,L)  are  computed  for  the  right 
side  of  the  subspan  and  0  is  incremented  for  the  next  subspan/point 
load. 

In  lines  149-165  renormalization  is  performed  as  described 
above,  involving  the  compilation  of  the  additional  1000  characters 
of  EMB21,  producing  a  total  length  of  about  5400  (about  5300  machine 
words).  If  further  modes  of  interest  remain  the  YES  branch  of  line 
177  is  taken  and  a  notice  printed  to  rerun  EMB15.  If  NM0DE=MINT(NMINT) 
=MSHP(NMSHP)  the  possibility  of  further  support  conditions  is  inves¬ 
tigated  in  statement  102;  the  NO  branch  signifies  problem  completion. 
The  YES  branch  causes  a  rerun  EMB15  notice  to  be  printed. 
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Assignment  of  Values  to  J,  K,  and 


EM'^25  4:54  SAT,  09/06/69. 

I  $FILE  EM8I0R,MB10'aE,MB10W5,MB15W1  ,M'?15W2/MBI5W3 

3  OIMENSI0N  \IPAR(7),PAR(7,17),IPAR<1,20),SUBSEG(17),SS(17), 

5  +  XN0SEG(17), C0EFM( 17),BMAT(23,4,4),STAMAT(24,4,1),C0RMAT(24,4, 1) 
7  +  VAL(6,70),DIST<70),A(4)^MINT(20) 

9  READ(1,^8)  ^SS, N0SUP, NMINT 

II  28  F0RMAT(4X,^3I5) 

13  READ(1,29)  (§S(  I )  ,  I  =1 . -^S) 

15  29  F0RMAT(4X,6F1O.5) 

17  READ(3,30)  (NPAR(I) , Ir 1 , 7) 

19  30  F0RMAT(17X,7I3) 

21  i^PARUNPARd) 

23  READ(3,31)  ( NPAR 1 , ( IPARC 1 , I ) , I = 1 , NPAR 1) ) ,SSUM 
25  31  F0RMAT(//7X,*I4//33X,E15.8) 

27  READ(2) 

29  D0  32  1=1 ,7 
31  MPRrMPARd) 

33  32  READ(2,33)  (PAR(I,  .1)  ,  J=I ,  NPR) 

35  33  F0RMAT(6X,4E15.8) 

37  READ(4,16)  MM0DE, MI NTCNMI NT) , 1 1, NATSP, ( C0RMAT( 1 , M, 1),M=1,4), 

39  +  (C0EFM( J), J=1 , NSS) 

41  16  F0RMAT(6X,I2,14X,I2,5X,I2,8X,I2/9X,4F4.1/(!1X,E15.8)) 

43  LIM=NSS+NPAR(1 ) 

45  REAn(5,18)  (STAMATC 1 ,M, H , «= I , t  I?  F0RMAT( /( 57X, E 1 4 . 8) ) 

47  D0  19  J  =  1  ,LIM 

49  REAr)(5,35)  A(l),A(2)j  35  F0RMAT<A6/A6) 

51  D0  19  K=1 ,4 

53  READ(5,20)  (BMAT<  J,  K,  L)  ,L  =  1  ,4)  ,  STAMAT(  J+ 1  ,  K,  1) 

55  20  F0RMAT(4E14.8,X,E14,8) 

57  19  C0NTINUE 

59  SEG=SSUM/20.J  LSUM=ITER:0 
61  D0  1  J=1 ,NSS 

63  IF(2-FIX<AINT(SS<J)/SEG)))  13,10,10 

65  10  XN0SEG<J)=2;  G0  T0  17;  13  XN0SEG( J) =AI NT(SS( J)/SEG) 

67  17  LSUMrLSUMf FIX<XN0SEG( J)) 

69  1  SUBSEG<J)=SS(J)/XN03EG(J) 

71  LIM2=LSUMf  NSS+2 

73  119  D0  50  K=l,4;  VALCK ,  1  )=STAMAT(  I ,  K,  1  ) ;  VAL(5, 1  )  =  VAL(3,  D* 

75  +PAR(7,1)/PAR(5,1);  50  VAL(6, 1 )=VAL(4, 1 )/PAR(6, 1 ) 

77  K=l;  L=2;  J=0;  DIST(1)=0. 

79  12  J=J+1;  DIST<L)=DIST(L-1) 

81  IF(J-IPAR(1, K))  3,4,3;  4  I F( NSS-IPAR( 1 , K) )  9,2,2 
83  2  D0  5  1=1,4;  5  VALCI ,L)=STAMAT( J+K, I , 1 ) 

85  VAL(5,L)  =  VAL(3,L)*PAR(7,  J)/PAR(5.J) 

87  VAL(6,L)=VAL(4,L)/PAR<5,J);  K=K+1;  G0  T0  8 
89  9  D0  15  1=1,4;  15  VALCI ,L) =STAMAT( J+K , I , 1) 

91  VAL(5,L)=VAL(5,L-1);  VALC6,L)=VAL(4,L)/PAR(6, J- 1 );  G0  T0  103 
93  3  D0  14  1=1,4;  14  VALCI ,L) =STAMAT(J+K-1 , I , 1 ) 

95  VALC5,L)=VALC3,L)*PARC7, J)/PARC5, J) 

97  VALC6',L)=VALC4,L)/PARC6,  J) 

99  8  SPANrO.;  LI  Ml  =X  N0SEGC  J)  ;  LSAVE=L;  M=0 

101  EI=PAPC4,  J)*PARC5,  J);  42  M=Mfl;  SPAN=SPANf  SU3SEQC  J) ;  L=L+1 
103  DISTCL)=D1STCLSAVE)+SPAN 

105  IFCLIMl-M)  44,40,44;  40  D0  43  1=1,4;  43  VALC I,L)=STAMATC J+K, I, 1 ) 
107  VALC5.L)=VALC3.L)*PARC7. J)/PARC5.J);  VALC6.L)=VALC4.L)/PARC6. J) 

Program  EMB25 
Figure  25 
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lOO  L=L+l;  G0  T0  12;  44  IF(PAR<3, J)-.0001 )  6,6,7 
111  6  VALCl  ,L)=VAL(  1  ,L-n  +  VAL(2,L-l)*SPAN 

113  +  +VAL(5,L-1  )*SPAM**2/(2.*EI)+VAL(4,L-1)*SPA!#*3/(6.*EI) 

115  VAL(2,L)=VAL(2,L-1)+VAL(3,L-1)*SPAN/EI+VAL(4,L-1)*SPAN**2/(2.*EI) 
117  VAL(3,L)  =  VAL(3,L-n  +  VAL(4,L-l  )*SPAN 

UR  VAL(4,L)=VAL(4,L-1);  VAL(5  ,L)  =  VAL(3 ,  L)*PAR(  7 ,  J)/PAR(5 ,  J) 

121  VAL(6,L):VAL(4,L)/PAR(6, J);  G0  T0  11 

123  7  TrC0EFMCJ);  ARG=T*SPAN 

125  U=(EXP(ARG)+EXP(-ARG))/2.;  V-C0S(ARG) 

127  X  =  (EXP(ARG)-EXP(-ARG))/2.;  Y=SIN(ARG) 

12R  W=U-V;  Z=X-Y 

131  VAL(1,L)=.5*(VAL(1,LSAVE)*(U+V)+VAL(2,LSAVE)*(X+Y)/T 
133  +  +VAL(3,LSAVE)*W/(EI*T**2)+VAL(4,LSAVE)*Z/(EI*T**3)  ) 

135  VAL(5,L)=.5*( VAL( 1 ,LSAVE)*W*EI*T**2+VAL(2,LSAVE)*Z*EI*T 
137  +  +yAL(3,LSAVE)*(U+V)  +  VAL(4,LSAVE)*(X+Y)/T)*PARC7,  J)/PAR(5,  J) 

139  VAL(6,L)=.5*(VAL( 1,LSAVE)*CX+Y)*EI*T*+3+VAL(2,LSAVE)*W*EI*T**2 
141  +  +VAL(3,LSAVE)*Z*T+VAL(4,LSAVE)*(U+V))/PAR(6,J) 

143  1 1  IF(LIMI-M)  41,41,42 
145  41  L=L+1 
147  G0  T0  12 

149  103  VALMAXzO.;  D0  104  L=1,LIM2 

151  104  VALMAX=AMAX1(ABS(VAL(1,L)),VALMAX) 

153  PRINT  177,VALMAX;  177  F0RMAT(7HVALMAX= , El  5,8) 

155  IF(ABS(ABS<VALMAX)-1.)-.001)  110,110,111 
157  111  RATI0=1  ./VALMAX 

159  G0  10(107,107,107,107,107,108,108,108,108,108),  NATSP 
161  107  M=l;  G0  T0  109;  108  M=4 

163  109  STAMAT(1,M,1)=RATI0*STAMAT(1,M,1);  ITER=0 
165  $USE  EMB21 

167  110  PRINT  178,  “REN0RMALIZED  M0DAL  VALUES  F0R  NM0DE=" , NM0DE 
169  178  F0RMAT(I2/7IHDIST  ALONG  BEAM  DEFLECTION  BENDING 

171  +  STRESS  SHEAR  STRESS) 

173  PRINT  179,  (D IST( L) , VAL ( 1 ,L) , VAL( 5 , L) , VAL(6 , L) , L= 1 , LIM2 ) 

175  179  F0RMAT<E15.S,4X,E15.8,4X,EI5.8,3X,E15.8) 

177  IF(MINT(NMINT)-NM0DE)  102,  102,  100 
179  100  PRINT, "RERUN  EMB15  FOR  NEXT  MODE" 

181  G0  T0  99 

183  102  IF(N0SUP-II)  98,98,101;  98  PRINT"  PROBLEM  COMPLETED  ;  G0  T0  99 
185  101  PRINT""RERUN  EMB15  FOR  NEXT  SUPPORT  CONDITION" 

187  99  CONTINUE 
189  S0PT  SIZE 
191  END 


LENGTH 

ABOUT  4400  CHARS. 
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Flaw  Chart  of  EMB25 


Figure  26 
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PRIMT  \ 

"PROBLEM  COMPLETED"  \ 


PRIIMT 

"RERUIM  EMB15  FOR  NEXT 
SOPPORT  COIMDITIOIM"  ■ 


109 


BTAMAT(1,M,1)= 

RATI0/STAMAT(1,M,1) 


ITER=0 


$USE  EMB21 


(DNCOO)^ 


VMmXz  .58490432E+01 
VaLMAX=  .10000119E+01 
RE^I0RMALIZED  M0DAL  VALUES 
DI3T  AL0NG  BEAM 


F0R  NM0DE= 


•OOOOOOOOE+Ol 
•OOOOOOOOE+Ol 
.17000000E+01 
.34000000E+01 
.34000000E+01 
. 42500000 E+01 
.51000000E+01 
.51000000E+01 
.61500000E+01 
.72000000E+01 
.72000000E+01 
.84000000E+01 
.96000000E+01 
.9S000000E+01 
.10300000E+02 
.1  1000000 E+02 
.1 lOOOOOOE+02 
.11450000E+02 
.1  1900000E+02 
.1 1900000E+02 
. 1 3400000 E+02 
.14900000E+02 
,14900000E+02 
.16600000E+02 
.18300000E+02 
.18300000E+02 
.18700000E+02 
.19100000 E+02 


DEFLECTION 
.17096813E+00 
.17096813E+00 
.13863570E+00 
.I0630383E+00 
.10630383E+00 
.901387S6E-0I 
.73981174E-01 
.73981 174E-01 
.54053567E-01 
.34258465E-01 
.34258465E-01 
.11963786E-01 
-.97661641 E-02 
-.97661  64  lE-02 
-.22083538E-01 
-.34076703E-01 
-.34076703E-01 
-.41601325E-01 
-.48  9  7  9  0  73E-01 
-.48979073E-0I 
-.72358416E-01 
-.93553 749E-01 
-.93553749E-01 
-.11602906E+00 
-.13816474E+00 
-.13816474E+00 
-.1433  1823  E+00 
-.14844937E+00 


,19100000E+02 
,20450000E+02 
.21800000E+02 
.21800000E+02 
.22850000E+02 
.23900000E+02 
,23900000E+02 
,24950000E+02 
.26000000E+02 
.26000000E+02 
,26550000E+02 
,27100000E+02 
,27100000E+02 
28650000E+02 
,30200000E+02 
,30200000E+02 
,31150000E+02 
,32100000E+02 
,32100000E+02 
32950000E+02 
33800000E+02 
33800000E+02 
34900000E+02 
36000000E+02 
3 60 00000 E+02 


-.14844937E+00 
-.16559389E+00 
-. 18245923E+00 
-.18245923E+00 
-.19537445E+00 
-.2081 1463E+00 
-.2081 1463  E+00 
-.22048921  E+00 
-.23235639  El-00 
-.23235639E+00 
-.21418935E+00 
-.15263101E+00 
-.15263101  EfOO 
.10843820E+00 
.39777000E+00 
.39777000E+00 
.567749  15  E+00 
.71891080E+00 
.71891080E+00 
.82807431 E+00 
.91054772E+00 
.91054772E+00 
.97764322E+00 
.10000119E+01 
.100001 19  E+01 


RERUN  EMB15  FOR  NEXT  MODE 


Eigenmode  and  Stmas 
Ualues  for  Second  Mode 


Figure  27, 

BENDING  STRESS 
.OOOOOOOOE+Ol 
.OOOOOOOOE+Ol 
.74800825E+00 
-.29724837E+00 
-.29724837E+00 
.40675727E+02 
22423  784E+01 
-.22423784E+01 
.53548799E+03 
.97672862E+03 
.97672862E+03 
.176521S8E+04 
.24719614E+04 
.24719614E+04 
.29774743E+04 
.34829734E+04 
.34829734E+04 
,38078715E+04 
.41328200E+04 
.41328200E+04 
.50961606E+04 
.60599847E+04 
.87235864E+03 
.  10284320  E+04 
.  1  1855985E+04 
.11855985E+04 
.12217342E+04 
.  125795 11 E+04 
.1257951 1 E+04 
. 134625 22E+04 
.14828485E+04 
.  14828485  E+04 
.13934587E+04 
,  13434725E+04 
.13434725E+04 
.  1  15534 23 E+04 
.  10107633E+04 
.84257230E+05 
.71825209E+05 
.60441613E+05 
.60441613E+05 
.  18985358E+05 
-.78795382E+04 
-.78795382E+04 
-.34404013E+05 
-.60931984E+05 
-.60931984E+05 
-,60954241E+05 
-.60984404E+05 
-.60984404E+05 
-.60991760E+05 
-.61003762E+05 
-.61003762E+05 


SHEAR  STRESS 
.OOOOOOOOE+Ol 
.OOOOOOOOE+Ol 
.42520182E-01 
.93335385E-01 
.93335385E-01 
.46033042E+01 
.97509838E+01 
.23696611  E+02 
.27327400E+02 
.31896800E+02 
.31896800E+02 
.33401524E+02 
.361 10421E+02 
.36110421E+02 
.36105009E+02 
.36103775E+02 
.36103775E+02 
.36095519E+02 
.360889 12E+02 
.32132053E+02 
.3 20 878 21 E+02 
.3206 1441  E+02 
.76947458E+01 
.76083735E+01 
.75397028E+01 
.75397028 E+01 
.75215105E+01 
.75039793E+01 
.  75^039  793E+01 
.35248485E+01 
-.44686274E+00 
-.5163  7540E+0.1 
-.90741855E+01 
-.12738841E+02 
-.12738841 E+02 
-. 17175322E+02 
-.21374701E+02 
-.35624502E+03 
-.39666836E+03 
-.42806465E+03 
-.42806465E+03 
-.44156728E+03 
-.46541 155E+03 
-.46541155E+03 
-.46526305E+03 
-.46516810E+05 
-.54045866E+00 
-.32738883 E+00 
-. 15306567E+00 
-. 15306567E+00 
-.68837025E-01 
.39500943E-02 
.39500943E-02 
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